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Abstract 

In this paper a new model method for describing of the electrostatic screening in single- 
component systems which is free of Debye-Hiickel's non-physical properties is presented. The 
method is appropriate for the determination of screening parameters in the case of the systems of 
higher non-ideality degree. The obtained screening characteristics are presented in a simple ana- 
lytic form. The presented results make basic elements of a method for determination of screening 
characteristics in dense two-component plasmas, which are discussed in Part2 and Part 3 of this 
work. 
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I. INTRODUCTION 



It is well known that, beside the existing strict methods, model methods are also used 
in all parts of plasma physics and in connection with most diverse problems. The purpose 
of every model method is to illuminate the real meaning of the problem considered in a 
physically acceptable, but much simpler way in comparison with the corresponding strict 



methods 
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28| . Probably the best known model methods is 



Debye-Hiickel's (DH method), which was developed for describing electrostatic screening in 
electrolytes (sf, but now for a long time it has been used in the plasma physics. Namely, 
a lot of electrostatic screening effects in plasmas are very often interpreted just in terms of 
such products of DH method as Debye-Hiickel's screening potential (DH potential), as well 
as Debye's radius and screening constant B, lo|, 12, 13, 16, 20, 27|. Because of that one 
could find these products of DH method in every course of plasma physics. 

When talking about DH method it is usual to include in this concept also the correspond- 
ing basic model of electrostatic screening in the considered system. According to [sj that 
model understands: 

• (al) the presence of an immobile probe particle, which represents one kind of charged 
particles in the real system (plasma, electrolyte); 

• (a2) treatment of free charged particles which fill the space around the probe particle 
as an ideal gas in the state of thermodynamical equilibrium; 

• (a3) treatment of the average summary electrostatic field as external with respect to 
that ideal gas. 

However, here we will adopt an agreement that in the further consideration "DH method" 
means exclusively the way of usage of the basic model, while this model itself will be con- 
sidered independently. Namely, although the basic model was introduced in 5|] together 
with DH method, this beneficent model has its own significance, since it can be used for 
development of another method for describing of electrostatic screening. In jsl the basic 
model was applied in a showy way which provided for DH method to possess the following 
very attractive features: 



(bl) - the procedure used is self-consistent, which means that all the existing quantities 
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are determined within the procedure itself and expressed through the basic plasma 
parameters (electron density N,,, temperature T, etc.), 

(b2) - the final results are presented by simple and compact analytical expressions. 

Certainly, it were these features that maintained the constant popularity of this method, 
alth oug h strict approaches based on classical or quantum statistical mechanics also exist 

However, the features (bl) and (b2) prevent from noticing the fact that DH method has at 
least two serious non-physical properties, masked by the applied mathematical procedure, 
which is discussed within this work. In this context one should remind to the problem 
connected with the screening constant for two-component systems. Namely, it was noticed 
that instead of DH constants, referring to the real two-component systems, the constants 
referring to the corresponding one-component systems were often used. Such a situation 
was probably first announced in {4], where DH screening constant of electron gas (on a 
positive charged background) was used instead of the one referring to the considered gase ous 



plasma. Later, gaseous plasma was treated in a similar way in several cases 
well as in the cases of very dense non-ideal plasmas 
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30|, as 



311] . This refiects the fact 



that one of mentioned non-physical properties manifests in every multi-component system 
independently of its non-ideality degree. It means that in the case of multi- component 
systems DH method in principle gives wrong results, even if the non-ideality degree is very 
small. Consequently, it was natural to expect that a new model method, free of all non- 
physical properties, but keeping positive features (bl) and (b2), would be more appropriate 
for plasmas of higher non-ideality degree, including dusty plasmas too. This was one of the 
main stimuli for this research whose purpose to develop a new model method which would 
be more adequate for description of inner-plasma electrostatic screening. 

An other stimulus was the fact that in the laboratory practice for a long time there 
has been a need for finding new screening characteristics of dense non-ideal plasma {Nf, = 
10^^ ^ lO^^cm"^, T = 1.5 ■ 10"^ ^ 3 ■ 10'^ i^) which would be different from Debye's radius 



and screening constant. It is illustrated by many papers 14 
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30|], where 



such new screening characteristics, introduced semi-empirically, were discussed. Besides, 
in some otherpapers such new characteristics were also used, although without any special 



discussion 
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3l| . Because of that, a need was evident for a new model method which 
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should be the generator of some screening characteristics appropriate for interpretation of 
experimental data. 

Already at the beginning of this research it was found that it is not possible to develop 
the searched model method by means of any corrections to DH one, and a completely new 
approach has to be found. However, we started from the analysis of the same basic model, 
keeping in mind that its exceptional properties (al)-(a3) allow including the inner-plasma 
electrostatic screening into consideration in the simplest way, and expecting that the model 
could offer some new possibilities. During this research it was found that DH method does 
not exhaust all the possibilities of the basic model, and that this model allows construction 
of a new method which posses positive characteristics (hi) and (h2), hut which is free of all 
the negative properties of DH one. The main aim of this work is just the presentation of 
this method, as well as several new screening characteristics obtained by its means. 

Due to the huge amount of problems which have to be studied in the general case of multi- 
component systems, only two simplest cases, i.e. the cases of single- and two-component 
systems, were considered in this work. The whole material has been arranged in three 
papers. The new model methods for single- and two-component systems are developed in 
Part 1 and Part 2, respectively, and the obtained new screening characteristic lengths and 
non-ideality parameters are presented in Part 3. The comparison of the results obtained by 
the developed method with the existing experimental data is also performed in Part 3. 

As the main objects of our researches we treated such two-component systems as fully 
ionized non- ideal gaseous plasmas (before of all hydrogen and helium ones), dusty plasmas, 
some electrolytes etc. However, in Part 1 we start just from single-component systems (an 
electron gas on the positive background etc.), since in that case the basic model itself is 
much simpler than in two-component case. Because of that, the procedure of elimination of 
one of the mentioned non-physical properties, which appears in all considered systems, can 
be developed in the simplest way. Consequently, in Part 2 it will be possible to concentrate 
whole attention to the elimination of other non-physical property of DH method. Also, 
we kept in mind the fact that a single-component system may always be treated as an 
approximation of a multi-component system. Therefore, the parameters of the corresponding 
single-component system could serve as an estimation in advance of the parameters of the 
observed multi-component system. 

The material presented in Part 1 is distributed in the next three Sections, as well as in the 
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three Appendixes. Section [TTl contains: description of the screening model, a critical analysis 
of DH method in the case of single-component system and stating the tasks precisely. In 
Section Hill the method developed for the single-component system is presented. Finally, 
Section IIVI contains results and discussion. 

II. THEORY ASSUMPTIONS 
A. Screening model 

A stationary homogeneous single- component system Sin is taken here as the initial model 
of some real physical objects. We will assume that: Sm is constituted by a gas of charged 
particles of only one kind and a non-structured charged background; the gas there is in 
an equilibrium state with temperature T and mean local particle density A^; the particles 
are point objects with charge Ze, where Z = ±1,±2, etc., and e is the modulus of the 
electron charge; the parameters Z, N and the background charge density pb satisfy the local 
quasi-neutrality condition 



Also, in the case of the electron gas it is assumed that the values of A^ and T allow its 
non-relativistic treatment. 

Electrostatic screening of a charged particle in the system is modeled with the help 
of the corresponding accessory system Sa which, accordingly to the basic model properties 
(al) and (a2), contains: a probe particle with the charge Ze, fixed in the origin of chosen 
reference frame (the point O), the charged background identical to that in Sm, and the gas 
of the free charged particles in the state of thermodynamical equilibrium with the same Z 
and T as in Sin- 

The gas of free particles in the system Sa will be characterized by: the mean local particle 
density n(r) = n(r; Z, T), the mean total charged density 



and the mean electrostatic potential $(r), where r = |r| and f is the radius- vector of 
the observation point. It is presumed that n(r) and p(r) satisfy the asymptotic boundary 



Pb + Ze-N = Q. 



(2.1) 



p{r) = pb + Ze - n{r) 



(2.2) 
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condition 



lim n(r) = N (2.3) 

1 — >oo 

and the condition of neutrality of the system Sa as a whole 

oo 

Ze + y p(r) ■ Ajrr^dr = 0. (2.4) 



Then, we will take into account that $(r) and p(r) are connected by Poisson's equation 

= -47r[Ze- (5(f) +p(r)], (2.5) 



where 6{f) is the three-dimensional (5-function 16|. It is presumed the satisfying of the 
boundary conditions 

lim $(r) = 0, (2.6) 

r— ►oo 

Iv^l < OO, (f = \im[^{r) - Ze/r]. (2.7) 
Since ip is the mean electrostatic potential in the point O, the quantity 

U = Ze-ifi (2.8) 

is the mean potential energy U of the probe particle. In an usual way we will treat U as an 
approximation of the mean potential energy of a free charged particle in the system Sin- 
In accordance with the basic model properties (a2) and (a3), as well as the boundary 
conditions fl2.3p and 02.61) . the condition of the keeping of thermodynamical equilibrium in 
the system Sa can be presented in the form 

fi{n{r),T) + Ze^r) = i2{N,T), (2.9) 

where fi{n, T) has a sense of the chemical potential of the ideal gas in the state of thermody- 
namical equilibrium with the particle density n and temperature T. Within this model just 
this equation is another one which, together with Eq. fl2.5p . provides the determination of 
charge density p(r) and the potential $(r). Let us emphasize that Eq. (I2.9|) i s ap plicable not 



only to the classical systems, but to the degenerated systems too (see [lO|, ll8|). However, 



one should keep in mind the difference between Eqs. ( 12. 5 p and (12. 9p . Namely, while Eq. ()2.5p 



is applicable in the whole space, Eq. (12. 9p is valid only in the area where usage of chemical 
potential fi{n{r),T) has a physical meaning. In the considered case, when the probe particle 
and free particles have the same charge, Eq. (12.91) is principally valid only in the region 



> 



3 



1/3 



(2.10) 



where is Wigner-Seitz's radius for the system Sm- Below will be treated as the probe 
particle self- sphere. 

In further considerations we will use the fact that equation ( 12.91) can be presented in the 
linearized form 



(2.11) 

n=N 



in the region of r, where 

\n(r) - N\ 

« 1- (2.12) 

Accordingly to the boundary condition Eq. (12.31) the part of space where the condition (I2.12p 
is satisfied always exists. 



B. The critical analysis of DH method 

The procedure of obtaining of DH solutions $D(r), pnir) and nci{r) for the electrostatic 
potential and the charge and particle densities in the single- component well as their 

properties, are described in Appendix |X1 The behavior of the reduced particle and charge 
densities, nD{r)/N and pnir)/ ph for one typical case is illustrated by Fig. [H This figure 
demonstrate apparent disadvantages of DH method: the negativity of the solution noir) in 
the region r < r^-"^ with the singularity in the point O; the existence of an additional non- 
physical region r^^^ < r^J. < r^, whose the sole role within DH method is to compensate the 
influence of the region r < rl~\ The distances r^^^ and r^~] are the roots of the equations 
( IA4p and (lASp . which always satisfy non-equalities < r|~^ < r^'^^ < r^. Apart of that, in 
Appendix |A] the attention is driven to the fact that for strongly non-ideal systems the direct 
manifestations of the non-physically behavior of no{r) appear also in the region r > r^. 

One can see that all non-physical properties of DH solutions are caused by the used 
procedure itself, which requires obtaining of the potential $D(r) as first. In this procedure 
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r [arbitrary units] 



FIG. 1: The reduced particle density n£,{r)/N and the reduced charge density pD{f')/pb in the 
case Kr<j = 1, where k is Debye's screening constant given by ()3.2p . 



$£)(r) is the solution of Helmholtz's equation (lAll) determined in the whole space, under 
the conditions (12.61) and (12.71) . since in the considered case the natural boundary conditions 
exits only r = oo and r = 0. In such a way it is completely neglected the fact that Eq. (lAip 
was obtained by means the equation (12. lip which is not applicable inside the probe particle 
self-sphere (0 < r < r^). 

C. What one should do in order to build a physically correct method? 

In accordance with above mentioned our main task will be the finding (within the basic 
model) of procedure which would provide satisfaction of the non-negativity condition 



n{r) > 0, 



< r < oo, 



(2.13) 
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as well as the application of the equation (12. lip only under the condition (12.120 . Also, we 
will keep in mind the fact that the basic equation (12.91) is correct only in the region (12.101) . 
wherefrom it follows that in the region r < the searched procedure mast not base on the any 
equation obtained from Ii2.9\) . including here the equation Ii2.11\) . It is clear that, contrary 
to DH procedure, the first aim of the searched procedure has to be determination of the free 
particle density n{r). 



III. THE METHOD PRESENTED 



A. The solution n(r) 

The region of large r: alternative procedure. In accordance with above request 
we will use Poisson's equation (12. 5p in order to obtain the potential <I>(r) for given charge 
density p(r). For this purpose we will use the expression (1C2p . In the form (]C2p the potential 
$(r) will be taken in the equation (12. lip . After the multiplication by Ze, it transforms in 
the integral equation of Volterra's type, namely 

CXJ 



p{r) = K' / p(r') -]r"dr', (3.1) 



where the screening constant k is given by the expression: 
1 



" (3 2) 

and is the one of the characteristic length which appear in considered model. Let us draw 
attention that k = in the classical case {dfi/dN = kT/N), and n = kt-f in the case of 
ultra degenerated electron gas (T = 0, dfi/dN = 2ep/3N), where 



1 

kd = — 
rD 



^ ' -Nl ,kt-f = =[—N] , (3.3) 



kT J rT-F \£f 



and and rjp-F are known Debye's and Thomas-Fermi radii, and Ef is the corresponding 
Fermi energy (see [lo|, Q). The principal significance of usage of just described procedure 
is caused by the fact that for any r in (13. ip only the region r < r' < oo appears. Because of 
that (13. ip can be treated within the region of r where the condition (I2.12p is satisfied. 
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The solution of the equation (13. ip will be found in the form: p(r) = S{r)/r. After that 
(13. ip gets the form 

oo 

S{r) = j Sir') (r' - r)dr'. (3.4) 

r 

In order to determine S{r) we will apply the operator d'^/dr'^ to both the left and the right 
side of (13. 4p and obtain the equation 

^ = .^5(r). (3.5) 

The general solution of this equation is S{r) = Aexp{—K,r) + Bexp{K,r). Since (13. 4p is 
satisfied only in the case B = 0, we have it that: p(r) = A ■ r^^exp(— Kr). On the base of 
this and (12. 2p we obtain, taking A = —Ze ■ a, the relation 

/ N exp(— Kr) 
n{r) = N-a-—^ (3.6) 

where is taken that a > 0. One can see that the alternative procedure, which is used here, 
provides analyzing of solution n(r) in the region where it satisfies the condition (I2.12p . The 
way of determination of the coefficient a in Eq. ( 13. 6p and obtaining n{r) in the whole region 
< r < cxD is described in details in Appendix |Bl 

The complete expression. From (13. 6p . (lB2p and (1B3P it follows that the complete 
solution n{r) is given by expression 

(AT AT / \ exp(— Kr) 
N - N -ro- exp(Kro) , r > ro, 
r (3.7) 
0, < r < ro, 

where the parameter ro can be determined from the condition (12. 4p taken in the form (13. 8p . 
In order to determine the radius Tq we will take into account that Eq. (12. 4p . after dividing 
with Ze, gives the equation 

oo 

J W - n{r)] ■ A-nr'^dr = 1, (3.8) 



which is especially discussed in Part 3. In the case when n(r) is given by Eq. (13. 7p . from 
(13. Sp it follows the equation 



;i + KTof = 1 + (Kr,)3, (3.9) 

10 



whose solution can be presented in two equivalent forms, namely 
'^o = '^s ■ ls{x), ro = r«;7«;(a;), 



(3.10) 



7s(x) = [(l + x3)l -l]/a;, 7«(a:) = (l + x=^)^ -1, (3.11) 
where the parameter x is defined by relations 

x = Krs = rs/r^. (3-12) 
From these expressions it follows that 

0<ro<r, (3.13) 
in the whole region < x < oo, and that 

limro = 0, lim Tq = r^. (3-14) 

Consequently, wigner-Seitz's radius represents the upper boundary for the radius tq. The 
physical meaning of the coefficients 7<j and 7^ is discussed in Part 3. The behavior of the 
solutions n(r) and p(r) is illustrated by Fig. [2] for urg = 1. 

B. The solutions p(r), $(r) and 
the potential energy U 



The expression for the solution p(r) is obtained by means of (12.11) . (12. 2p and (13.71) . and 
it is presented here in the form 

7,r /X exp(— Kr) 

I -ZeN ■ro-exp{Kro) , r > Tq, 

p(r) = { r (3.15) 

-ZeN, < r < ro, 

where tq is given by Eqs. (I3.10p and (13.121) . The corresponding expression for $(r), obtains 
by means of (lC2p and (IC4p in the form 



Ze 



exp(— /tr) • x(2;), r > tq. 



r 
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FIG. 2: The reduced particle density n{r)/N and reduced charge density p{r)/p{y in the case 



Kr.o = 1. 



where the factor x{^) is given by expression 



■ exp 



(l + 



(3.17) 



The behavior of xi^) is illustrated by Fig.[3]and discussed in Appendix[Bl The point x = 7^ 
is shown in this figure since in this point ro = r^. 

Accordingly to Eq. (12. 8p determination of the potential energy U of the probe particle 
requires knowledge of the potential ip defined by Eq. (12.71) . By means of Eqs. (13.15l) - (l3.12p 
and Eq. flC3l) one obtains that 

From here and Eqs. (13.111) it follows that the potential energy U can be presented in two 



12 




FIG. 3: The behavior of the function x{^)i defined by (j3.17p . In the point x = the equaUty 
ro = is vahd. 



equivalent forms 



3 (l+a;=^)^ -1 ^ {Zef 
2 ^ ^ 



t/ = t/.-- ' ^ , U^=-'—^, (3.19) 



t;.„,.il±4l^, U..-l-^, (3.20) 



where f/^ is DH potential energy of the probe particle (see Appendix |Al), and Us is the 
potential energy of this probe particle only in the field of the charged background which 
fills its self-sphere (0 < r < r^). The parameters f/^ and Us represent the corresponding 
boundary potential energies, since from (13.191) and (I3.20p it follows that 

lim U/U^ = 1, lim U/Us = 1. (3.21) 

X— >0 X— >oo 

The behavior of the the ratios U/U^ and U/Ug is illustrated by Fig. HI 
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FIG. 4: The behavior of the ratios U /Uf^ and U /Ug as functions of x. The point x = 3/2 represent 
an arbitrary boundary between the regions of weak and strong non-ideahty. 

IV. RESULTS AND DISCUSSION 

The expressions fl3.7p - fl3.12p and fl3.16l) - fl3.20p show that the obtained solutions ?7,(r), 
p(r) and $(r) satisfy all conditions introduced in Section [Tll can be applicable to the single- 
component systems for any /tr^ > 0. It is important that from Appendix [B] it follows that 
within the basic model (al)-(a3) the used alternative procedure is unique one which provides 
that the solution n{r) satisfies both the neutrality condition (12. 40 and the non-negativity 
condition f l2.13p and and posses features (bl) and (b2). 

Accordingly to Eqs. (13.70 and (13.100 - 03.120 . we have always that n(r = 0) = 0. In the 
case of the classical system such a condition has to be satisfied from energetically reasons 
for any Z and T. However, in the case of degenerated system, when the free particles have 
to be treated as quantum mechanical objects, n(r = 0) is the small positive quantity. In 
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this case the mentioned condition practically solve the problem of the elimination of DH 
negative singularity, since provides the transition from ^^(r = 0) = — oo to n{r = 0) = 0. 

From the expressions fl3.17p . (I3.16P and (13.191) one can see that the values of the potential 
$(r > Vg) and the probe particle potential energy U are closed here to the corresponding 
DH values > fs) and in the region urg <^ 1. This fact reflects the specificity of 

the single-component systems where the screening constant k is the same in both DH and 
presented methods. However, the behavior of $(r) and U essentially differs from the behavior 
of $j:)(r) and in the single-component case for large values of nvg. So, accordingly to 



Just such a case should realize in dusty plasmas containing dusty particles with the charges 
Ze, where \Ze\ ^ 1 and the temperatures T ^ lOOOi^. The figure H] suggests that for U it 
is suitable to use Eq. f l3.19p for x < 3/2, and Eq. f l3.20p for x > 3/2. The point x = 3/2 can 
be interpreted as an arbitrary border between the region of weak non-ideality (DH region) 
and the region of strong non-ideality (non-DH region). 

One of the main results of the method which is developed in this paper is appearing 
of three new parameters tq, 7^ and 7^ which are given by Eqs. (I3.10p . (13. lip and (13.121) . 
The parameter ro represents the radius of the sphere centered in the probe particle which 
is classically forbidden for the particles from their environment, while js and 7^ can be 
interpreted as some kind of non-ideality parameters. Full sense of these quantities will be 
discussed in Part 3. Also, in Part 3 the results of this paper will be used for obtaining 
of other relevant screening parameters. Finally, let us draw attention that the alternative 
procedure, described in this paper, will be applied in Part 2 of this work in connection with 
two-component systems. 
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Eq. (Km and (ICTIl . U ^ Us when nr^ > 1 and 




(4.1) 
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APPENDIX A: DH SOLUTIONS 



1. The procedure and expressions. 

In DH procedure the free particle density is expressed by means of equation (12. lip over 
the electrostatic potential and in such a form used in Poisson's equation (12.51) . Then, by 
means of (12.11) and (12. 2p . Eq. (12. 5p transforms to Helmholtz's equation: 

V^$(r) = fi:^<l>(r), (Al) 

which applies in the whole space < r < oo, neglecting the conditions (12.101) and (12.121) . 
The solution $z)(r) is determined by the boundary conditions (12.61) and (12.71) . After that, 
DH charge and particle densities poif) and nD{r) are obtained in the whole region r > 0, 
by means of ( 12.21) and ( 12. lip . These solutions are given by 

^oi''") = Ze ■ exp(— Kr), 

Pz)(^) = —{ZeK^/ATx) • r^^ exp(— /tr), (A2) 
noir) = N — [k^ /Atx) ■ exp(— /tr), 

where n is given by (13.21) or (13. 3p . DH values ipn and Un of the potential in the point O and 
the probe particle potential energy are determined from (12.70 . (12.80 and (]A2p and given by 

LpD = -Ze-K, Ud = -{Zef ■ K = -{Zef /r,,, (A3) 

which are very often used in plasma physics. One can see that DH value Ur, = U^, where 
f/^ denotes the boundary value of the potential energy in the expression (I3.19p . 



2. The solution n£i(r): the region r < r^. 

From the expression (]A3P it follows that nD{r) < in the region r < r^^ \ where r^^^ is 
the root of the equation 

nnir) = 0, (A4) 

which exists for any k > 0. Then, the same expression for n£,{r) shows that the existing of 
the region where n£)(r) < is compensated within a wider region r < r^J, where r^^^ is the 
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root of the equation 



r 



I 



nr,{r') ■ Anr'^dr' = 0, 



(A5) 





which also exists for any k > 0. From the neutrahty condition (12.41) it follows that always 
^out < ^s- All mentioned is illustrated by Fig. [TJ 

3. The solution n/)(r): the region r > r^. 

We will take into account that in accordance with flA2l) and (13.121) the condition fl2.12p 
with n{r) = n£,{r) can be presented in the form 



From here it follows that the quantity h£){r) monotonously increases when f decreases, for 
each fixed x. Therefore, in the region r > r^, the quantity h£){r) reaches its maxima at 
r = Tg, which it means that it is enough to consider the behavior of the quantity hni'^s)-, 
which is given by 



as a function of x. From this expression one can see that hu{rs) = at x = and x = oo, 
as well as that this quantity reaches its maxima at x = 2. Since 



we can consider that n£,{r) satisfies the condition (I2.12p in the region r > r^, for any Kr^ > 0. 
However, the decreasing h£){rs), when x increases in the region x > 2, means that in this 
region the behavior of n£){r) becomes non-physical in the whole space. 

APPENDIX B: THE BEHAVIOR OF THE SOLUTION n(r) 
1. The region r < r^: the extrapolation procedure. 




(A6) 




(A7) 




(A8) 



Because of the discussion in lll Cl in the region r < is not possible using of any procedure 
based on the equation (12. 9p . we have to return for a moment to the solution n£)(r) and 
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consider it from the aspect of the possibihty to obtain it avoiding the procedure described 
in lA 1[ Namely, noir) can be treated as a result of an extrapolation of the expression 
(13. 6p in as wide as possible area (0 < r < oo) allowed by the condition fl2.4p . Such a way 
of extrapolation provides that the coefficient a in fl3.6l) takes just DH value h'^/Att. This 
procedure is consistent one since rinir), accordingly to Eq. (lASp . satisfies the condition 
f l2.12p in the region r > for any nrg > 0. 

The solutions obtained in a similar manner are already known in physics. For example it 
is enough to mention Slater's and Bates-Damgaard's wave functions which have successfully 
been used in atomic physics , l25|, |26|]. The applicability of that solutions was caused 



by adequately chosen boundary conditions. However, in the case of solution ^^(r) such 
conditions were not used, what caused its negativity in the region < r < r^^\ 

In connection with this one should keep in mind that noir), independently from its non- 
physicality in the region < r < r^~], where r^~^ < r^'j < Vg, provides acceptabihty of the 
potential $£i(r), as well as the screening characteristic length r^, for weakly and moderately 
non-ideal single-component systems {kTs < 1). This fact suggests that an adequate proce- 
dure of extrapolation of the expression ( 13.60 . which excludes appearing of negative values of 
the solution n{r), as well as its non-physical behavior in the point O, would provide the ap- 
plicability of n(r) not only to weakly and moderately non-ideal systems, but to the systems 
of higher non-ideality {nvg ^1). 

Accordingly to above mentioned the procedure of obtaining of n{r) has to be continued 
by the extrapolation of the expression (13. 6p in as wide as possible region of r, namely 
ro < r < oo, allowed by the condition (12.130 . From ( 13. 6p it follows then that ro is the root 
of equation 

N -a- exp(-Kr) /r = 0, (Bl) 
and consequently 

a = iV ■ ro exp(/tro). (B2) 

In the region < r < tq the solution n{r) has to be continued by means of equality 

n{r) = 0, < r < ro, (B3) 

which provides that this solution at least does not increases when r ^ 0. The corresponding 
form of the complete solution n{r) is given by the expression (13. 7p . The parameter ro in this 
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expression has to be obtained from the neutrahty condition in the form (13. 8p . The same 
condition shows that it is always ro < r^. 

The complete procedure of obtaining of n(r) is consistent one since n(r), given by (13 .Tp 
and fl3.10p - fl3.12p . satisfies the condition (12.120 in the region r > for any nr^ > 0, as 
it follows from flBSp . Finally, it is very important to draw one's attention that described 
extrapolation procedure is unique one which provides that the solution n(r) is self-consistent 
and simultaneously satisfies all conditions from Section [Til including the additional condition 
fl2.13p . Namely, as it can be shown, any other extrapolation procedure causes the appearance 
of at least one parameter that cannot be determined within that procedure itself. 



2. The region r > r^. 

From (12. ip . (13. 7p . fl3.10p and f1A2p it follows that the product Nroexp^uro) in the expres- 
sion for n{r) differs from the coefficient k^/Att in the expression for noir) only by the factor 
x{x). The behavior of x{^): given by fl3.12p and fl3.17p . is presented in Fig. [3l One can 
see that in the region < x < 2 this function increases from 1.0 to a value that is close to 
1.2. From this it follows that for urg < 2 the solution n{r) has to automatically satisfy the 
condition fl2.12p . In the region urg > 2, where the non-physical properties of n£){r) appear, 
we have to directly analyze the left side of Eq. f l2.12p . From here and (13. 7p it follows that 

h{rs) = {x /3) ■ exp(-x) ■ x{x), h{r-s) = , (B4) 

and, for the difference of holrs), monotonously increases whit the increasing of x in the 
whole region x > 0. Then, on the base of fl3.17p and f1B4p . we obtained that 

lim h{rs) = e'^ . (B5) 

X— >oo 

Above mentioned means that n{r) is free of the non-physical properties, which have nD[r), 
and satisfies the condition f l2.12p in the region r > Ts for any kTs > 0. 
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APPENDIX C: THE SOLUTION ^>(r) AND THE POTENTIAL 99 



As it is known, the solution of Poisson's equation (12.51) . which satisfies the boundary 
conditions (12. 6p and (12.71) . is given by the expression 



<l>(r) 



r 00 



(CI) 



This expression, by means of the electro-neutrahty condition (12.41) . gets the form 



<l>(r) = — 

r 



+ 



p(r')47rr'^(ir' 
wherefrom it follows the expression 
$(r) = -471 J p{r') - ^r'^dr'. 



Pir'] 



A-nr''^dr' 



(C2) 



That is one out of the two expressions for the potential $(r) which are used within the frame 
of this work. The basic feature of (1020 is that for each r > the potential $(r) is expressed 
only by p(r') from the region r' > r. 

From the definition (12. 7p of the potential <f and (]C1|) it follows the expression 



p{r) 



Anr^dr, 



(C3) 



which uses in this work for the determination of (f. 

Finally, by means of (ICip and (1C3P one obtains another expression for the potential $(r), 
namely 



$(r) 



Ze 



r 



(C4) 



The basic feature of ( lC4p is that for each different r = ri andr = r2, the difference 
[$(r2) — ^'(?"i)] is expressed only by p(r') in the region < r' < max(ri,r2). This fact 
will play a significant role in Part 2. 
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